Identification of Bursting Water Maser Features in Orion KL 
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ABSTRACT 

In February 2011, a burst event of the H2O maser in Orion KL (Kleinmann- 
Low object) has started after 13-year silence. This is the third time to detect 
such phenomena in Orion KL, followed by those in 1979-1985 and 1998. We 
have carried out astrometric observations of the bursting H2O maser features 
in Orion KL with VERA (VLBI Exploration of Radio Astrometry), a Japanese 
VLBI network dedicated for astrometry. The total flux of the bursting feature at 
the LSR velocity of 7.58 km s~^ reaches 4.4x10^ Jy in March 2011. The intensity 
of the bursting feature is three orders of magnitudes larger than that of the same 
velocity feature in the quiescent phase in 2006. Two months later, another new 
feature appears at the LSR velocity of 6.95 km s~^ in May 2011, separated by 
12 mas north of the 7.58 km s~^ feature. Thus, the current burst occurs at 
two spatially different features. The bursting masers are elongated along the 
northwest-southeast direction as reported in the previous burst in 1998. We 
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determine the absolute positions of the bursting features for the first time ever 
with a submilh-arcsecond (mas) accuracy. Their positions are coincident with the 
shocked molecular gas called the Orion Compact Ridge. We tentatively detect 
the absolute proper motions of the bursting features toward southwest direction. 
It is most likely that the outflow from the radio source I or another young stellar 
object interacting with the Compact Ridge is a possible origin of the H2O maser 
burst. 

Subject headings: ISM: individual objects (Orion KL) — ISM: molecules — 
masers — radio lines: ISM 



Introduction 



In the context of astronomical masers, an enormous outburst of the II2O maser at 
22 GHz (A=1.3 cm) in the Orion KL (Kleinmann-Low object) region is one of the most 
enigmatic phenomena. The first II2O maser burst in Orion KL was discovered in 1979 and 



the active phase continued until 1985 with several fiare-up events (IMatveenko et al.lll988 



Garay et al.lll989l ). Based on the VLBI (very long baseline interferometer) observations, the 
bursting maser fea ture was revealed to hav e east-west elongation which was interpreted as a 
circumstellar disk (IMatveenko et al.lll988l ). Followed by the active phase for about 6 years, 
the quiescent phase lasted from 1985 to 1997. In t his phase, a jet-like mas er feature with 
the northwest-southeast elongation was identified (IMatveyenko et al.Ml998l ). The second 
burst was detected at the end of 1997 and monitoring observations were carried ou t with 
single-dish radio telescopes and VLBI (jOmodaka et al.lll999l : IShimoikura et al.ll2005l ). The 
structure in the second burst phase was consistent with that in the adj acent quiescent phase 
(jOmodaka et al.lll999l : IMatveyenko et al.ll2004j : IShimoikura et al.l 120051 ). The maximum flux 
density flared up to an order of 10^ Jy in both burst phases. The line-of-sight velocity with 
respect to the local standard of rest (LSR) of the bursting maser was ~ 8 km s~^ and was 
unchanged in both phases, suggesting a common origin. 

Two different scenarios hav e been proposed for this outburst; one is relate d to the exis- 



tence of a circumstellar disk (e.g. Matveenko et al. 19881 : Shimoikura et al. 2005 ) and another 



is the jet /outflow driven by y oung stellar object (YSO) in the Orion Nebula (iGaray et al. 



19891 : IMatveyenko et al.ll2004l ). However, no evidence of either circumstellar disk or the 



jet/outflow system associated with the bursting maser has been identified, and hence, the 
origin of this maser burst is still unclear. 



The third-time burst of the H2O maser in Orion KL has started since February 2011 
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( iTolmachevI 1201 ll ). The LSR velocity of the bursting maser is again 7.8 km s~^. The burst 
appears to be possibly periodic with an interval of about 13 years (1985, 1998, and 2011), 
which could support the hypothesis of the common origin. 

In order to explore the mechanism of this maser burst, we have started astrometric 
observations of the bursting H2O maser in Orion KL with VERA (VLBI Exploration of 
Radio Astrometry), a Japanese VLBI network developed for astrometry. Here we report the 
result of our absolute position measurement of the bursting maser features with a submilli- 
arcsecond (mas) accuracy. It allows us to measure the proper motions of the bursting 
features. Further detailed studies on the structure and time variation of the bursting maser 
features will be undertaken. 



Observations 



Observations of the H2O maser (6i6-523, 22235.080 MHz) in Orion KL were carried out 
on 2011 March 09, May 01, and May 17 with VERA. The weather conditions were good 
at all of the four stations (see Figure 1 of iPetrov et al.l 120071 ). and the typical system noise 
temperature was 100-200 K. The maximum baseline length was 2270 km and the uniform- 
weighted synthesized beam size (FWMH) was 1.7 masxO.9 mas with a position angle of 
143 degrees on average. We employed t he dual beam ob servation mode, in which Orion KL 
and an ICRF source J054138. 0-054149 (IMa et al.lll998l ) were observed simultaneously. The 
data were recorded onto magnetic tapes at a rate of 1024 Mbps, providing a total band- 
width of 256 MHz in which one IF channel and the rest of 15 IF channels with 16 MHz 
bandwidth each were assigned to Orion KL and ICRF J054138. 0-054149, respectively. For 
the H2O maser lines, a spectral resolution was set to be 15.625 kHz, corresponding to a 
velocity resolution of 0.21 km s~^. A bright continuum source, ICRF J053056. 4+133155, 
was observed every 80 minutes for bandpass and delay calibration. Calibration and imaging 
were performed using the NRAO Astronomical Image Processing System (AIPS) software 
package. 

In addition to the present observ ations, we reanaly zed the archive data of the previous 
monitoring observations of Orion KL (IHirota et al.ll2007l ) conducted on 2006 July 15 in order 
to compare the results of the new ob servations. Details of the observation and data analysis 
are described in iHirota et al.l (120071 ) . 
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3. Results and Discussion 

3.1. Images and spectra of the bursting maser features 

The H2O maser burst is detected at the peak LSR velocity channel of 7.58 km s~^ in 
the first epoch of observation on 2011 March 09, as shown in Figure [H The total flux density 
is (4.4±0.3) xlO^ Jy. The peak velocity is apparently shifted toward lower velocity ch annel 



compared with that in the previous burst event, 7.64 km s~^ (jShimoikura et al.ll2005l ). The 
correlated flux, (2.2 ±0.2) x 10^ Jy, recovers 50% of the total flux density. Due to the lack of 
short baselines in VERA, extended emission component would be resolved out. Two months 
later from the flrst epoch, another new velocity component at 6.95 km s~^ appears in 2011 
May 01 while the flux density of the 7.58 km s~^ component gradually decreases. The total 
flux density of the 6.95 km s~^ component becomes comparable with that of the 7.58 km s~^ 
component on 2011 May 17. 

We have successfully obtained phase- referenced images of the bursting H2O maser fea- 
tures. The position, peak intensity, and size of each feature are listed in Table |T1 Note that 



due to the loss of coherence caused by the atmospheric zenith delay residual (IHirota et al. 



20071 . and references therein) in the phase-referencing calibration, the peak flux density is 
degraded by a factor of 0.5-0.9 compared with the results obtained from the conventional 
self-calibration method (columns 6 and 7 in Table [1]). Figure [2] shows the phase-referenced 
images of the bursting maser features at the LSR velocity channels of 6.95 km s~^ and 
7.58 km s~^. In the flrst epoch in March 2011, the 7.58 km s~^ feature shows single-peaked 
structure while it splits into double peaks in May 2011. In addition, another spatially dis- 
tinct feature at the LSR velocity of 6.95 km s~^ appears at 12 mas north as also seen in the 
spectra (Figured]). Therefore, the current maser burst is turned out to occur at two spatially 
different features. 

The maser features elongate along the northwest-southeast direction. Results of the 
two-dimensional Gaussian fltting on the phase-referenced images are also summarized in 



(Hirota et al. 


2007; 


Kim et al. 


2008) 



of the bursting maser features are sometimes resolved along their major axes (i.e. those 
with the major axes of larger than 2 mas). The 7.58 km s~^ feature shows a double-peaked 
structure with a separation angle of 2.7 mas on 2011 May 01 and 17 as shown in Figures 
Mjo) and (c). The size of the maser feature is comparable to or slightly smaller than that 



found in the previous burst in 1998 (1.8 masxO.8 mas-5.3 masxl.5 mas; IShimoikura et al. 
2OO5I ) . The brightness temperature reaches 2.3x10^^ K. The velocity width of the total power 



spectrum of the bursting maser is 0.76±0.07 km s ^ (FWHM) for the first epoch. When 
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we fit the velocity width of each feature at different velocity and position derived from the 
Gaussian fitting of the self-calibrated images (Tabled]), we obtain the average velocity width 
of 0.61±0.03 km s~^. These velocity wid ths are broader than that in the previous burst, 
0.40-0.48 km s'^ dshimoikura et al.lboosi ). 



3.2. Astrometry of the bursting maser features 



In the previous active phases, high- resolution observations were conducted with the VLA 



and VLB! ( 


Matveenko et al. 


1988 


Garav et al. 


1989: 


1999; 


Matvevenko et al. 


2004 


: Shimoikura et al. 


2005) 



Greenhill et al. 1998; Omodaka et al. 



racy was ~50-100 mas even in the most extended A- configuration (e.g. lGreenhill et al.lll998l ). 
Although VLBI observations achieved much higher relative position accuracy of 0.05-0.1 mas, 
the absolute position information was usually lost by the self-calibration analysis. 

In contrast, we have carried out phase-referencing VLBI astrometry, in which one can 
measure the position of the target maser source with respect to the extragalactic position 
reference source. For instance, the accuracy of absolute position measurement with VERA 
is reported to be 0.36 mas and 0.74 mas in right ascension and declination directi ons, respec- 
tively, for the previous H2O maser observations in Orion KL (IHirota et al.l 120071 ). The main 
source of the astrometric errors in the VERA observations is the variation of the structures 
of maser features as well as the atmospheric zenit h delay residual due to the troposphere 
and the positional accuracy of the reference source (IMa et al.lll998l ). 



The observed maser positions are summarized in Table [U This is the first time ever to 
measure the absolute positions of the bursting maser features in Orion KL with the submilli- 
arcsecond accuracy. The current position of the bursting maser is shifted by +70 mas and 
-200 mas in right a scension and declinat ion, respectively, with respect to that of the previous 
burst in 1983 (e.g. iGreenhill et al.lll998[ ). probably because of much worse astrometric accu- 
racy of previous observations than that of our pre sent result. On the oth er hand, the position 
of the bursting maser feature measured in 1998 (jOmodaka et al.lll999[ ) is different by 2" in 
declination. This is due to the lower positional accuracy of an order of ~1" in the fringe- 
rate mapping method in particular for declination which is typical for the interferometric 
observations of target sources close to the equator such as Orion KL. 

We have measured the absolute proper motions of three bursting maser features as 
shown in Ta ble [2] and Figure |3l Taking into account the annual parallax of Orion KL 
of 2.39 mas (IKim et al.l |2008[ ). we obtain the absolute proper motions of ~4-9 mas yr~^ 
or 8-18 km s^^ toward southwest direction. At first glance, they seem inconsistent with 
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the apparent position movement in the phase-referenced images due to the relatively large 
parallax approximately toward the eastward direction. Although the present results contain 
large uncertainties due to the short monitoring period (~3 months), the absolute proper 
motions could be attributed to the outflow in the Orion KL region as discussed below. 

We tentatively identify the possible counterpart of the bursting maser in the quiescent 
phase data on 2006 July 15 at the peak LSR velocity channel of 7.37 km s~^, according to the 
similarity in the line-of-sight velocity, position (within 100 mas from that in March 2011) and 
the structure (NW-SE elongation). The position is also listed in Table [TJ The peak intensity 
of the phase-referenced image is only 19.6±1.4 Jy beam~^ on 2006 July 15 (see Tabled]). 
The peak intensity derived from the conventional self-calibrated map is 39.6±0.5 Jy beam~^. 
The intensity of the bursting feature is three orders of magnitude greater than that of this 
feature in the quiescent phase. 

Judging from the similarities in peak LSR velocity, position, and structure among the 
detected features in 2011, 2006, and earlier, they could be attributed to the common origin. 
However, it is unlikely that the bulk of the bursting gas clump would be identical throughout 
every burst and quiescent phases because the time span of observations (>5 yrs) is longer 



than the typical life time of the H2O masers in Orion KL (~1 yr; see Figure 1 in iHirota et al. 



20071 ). In fact, we can see significant time variation of velocity and spatial structure of the 
bursting maser features even within the present monitoring observations for only two months. 
If the gas clump of the bursting maser feature would be identical to that of the possible 
counterpart detected in 2006, one of the features should show an absolute proper motion of 
(-1-8 mas yr^^, -17 mas yr~^), which corresponds to the position movement by (-1-34 mas, 
-80 mas) in right ascension and declination, respectively, during the 5-years separation. We 
did not detect such a large proper motion in our astrometric data from March to May 2011. 

We also find two other maser features at the same velocity range ~8 km s~^ in the 
quiescent phase in 2006. These positions are shifted by (-1.0",-3-l") and (-1-1.7", -5.4") from 
the bursting maser features. If either of them was identical to the bursting feature, the 
proper motion would be unrealistically large, 1400-2400 km s~^, according to the position 
difference during the 5-year separation in our VERA observations. Therefore, we conclude 
that no maser feature detected in 2006 is identical to those in the bursting phase in 2011. 



3.3. Origin of the bursting maser features 



In this section, we discuss a possible powering source of the bursting maser features. 
One of the most plausible interpretations is related to the outflow from YSO(s) in the 
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Orion KL region as proposed by iGaray et al.l (119891 ) . It has long been established that 
the H-;.0 masers in Ori on KL are excited by interaction with outflow and ambient dense 
gas (IGenzel et al.l Il98ll ) . The most possible p owering source is the radio source I which 
drives northeast-southwest low- velocity outflow (jPlambeck et al.ll2009l ). The bursting maser 
features are located at 8" or 3400 AU southwest from source I, and is coincident with the 
interacting region b etween the outfl ow and a dense ambient gas, which is called the Orion 
Compact Ridge (e.g. lLiu et al.ll2002l ). Although the elongation of the bursting maser features 
(PA=130-150 degrees) are perpendicular to the direction to source I, the maser features could 
trace the shocked layer between the low-velocity outflow and the Compact Ridge. The H2O 
maser burst may occur at the different part of the shocked layer when the episodic or possibly 
13 year-periodic outflow from source I interacts with the Compact Ridge. The magnitude 
of the proper motions of the bursting maser features are 16-20 km s~^ pointing toward 
the west to southwest direction when we subtra ct the absolute pro per motion of source I 
(6.3 mas yr~^, -4.2 mas yr~^) recently reported by lGoddi et al.l (120 1 ll) . Interestingly , they are 
roughly consistent with the low- velocity outflow from source I (IGenzel et al.lll98ll ) although 
this should be confirmed by the further VLBI monitoring observations. 

Nevertheless, it is still unclear why only the 8 km s~^ features show such anomalous out- 
burst. It should be noted that the current burst event cannot be explained by the local phe- 



nomenon around each feature with the scale of a few milliarcseconds (e.g. iMatveenko et al. 



19881 : IShimoikura et al.ll2005l ) because the maser burst is occurring at two distinct features 
separated by 12 mas and at different velocities. One of the possibilities is that powerful 
outflows from source I may interact with the dense core associat ed with a pre-existing YSO 
in the Compact Ridge region as suggested by I Gar ay et al.l (119891 ). What makes the bursting 
features al ways at about 8 km s~^ would be that it is the rest velocity of the Compact 
Ridge (e.g. iLiu et al.ll2002l ). The low- velocity outflow hitting at any fortuitous place in the 
Compact Ridge may amplify the maser emission most efficiently at the rest velocity, leading 
to such features. According to the SIMBAD database, there are number of closer infrared, 
radio and X-ray sources around the bursting maser features than source I. The nearest YSO 
candidate is the radio continuum source R, which is 0.5" southwest from the bursting maser 
features . A millimeter an d infrared counterpart, HC438, is also associated with the same 
source ( Eisner et al. 20081). Submil l imeter continuum emission peaks are detected within 



their beam widths (ITang et al.ll2010l : IZapata et al.l 201 ll) . Although the mass of the submil- 
limeter core is estimated to be ~4 M0 (ITang et al.ll2010l ). the nature of this source, whether 
it is a self-luminous YSO or not, is still unclear. Thus, we cannot definitely identify one of 
these sources as a powering source of the bursting maser features only based on its proximity 
at this moment. 
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4. Summciry 

We have successfully measured the absolute positions of the bursting H2O maser features 
in Orion KL with the VLBI astrometry technique using VERA. This is the first time ever 
to identify the bursting maser features with better than the milli-arcsecond accuracy. The 
burst event is found to be occurring at two spatially different maser features. These features 
are located in the interacting region of the outflow from the radio source I and the Orion 
Compact Ridge. Although the spatial and velocity structure of the bursting masers are 
variable with time, the position and structure of the bursting masers are similar to those 
in the previous burst phases, implying a common origin. We suggest that the outflow from 
source I or another YSO in the Compact Ridge region could be a possible powering source 
of the bursting maser features. 

If the current burst phase persists for one year as in the case of previous bursts, we 
will be able to determine the proper motions of the bursting maser features more accurately 
through our VLBI astrometry with VERA. They will be crucial to reveal dynamical structure 
of the masing region. For instance, we will provide novel information on the VLBI images 
and proper motions to construct three-dimensional dynamical model of the outflow. If the 
bursting masers are associated with a circumstellar disk, we will be able to derive the mass 
of the central protostar by directly measuring the rotating motion. In addition, high spatial 
resolution observations with ALMA and EVLA of the thermal molecular lines and continuum 
emissions from centimeter to submillimeter wavelengths are still indispensable to identify the 
powering source of the bursting maser along with the velocity structure of possible outflow 
from it. All of these results are essential to shed light on the emission mechanism and origin 
of the burst of the maser features in Orion KL, as well as the driving source(s) of complex 
outflows in this region. 

We thank the staff of the VERA stations to support our project, Shuji Deguchi, Satoshi 
Yamamoto, and anonymous referees for useful comments. This work is supported in part 
by The Graduate University for Advanced Studies (Sokendai), and has made use of the 
SIMBAD database, operated at CDS, Strasbourg, France. 
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Table 1. Properties of the bursting maser features 
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''Feature IDs are shown in Figure [2] 

''Peak intensity derived from the phase-referencing method. 
'^Peak intensity derived from the self-caUbration method. 



Note. — Reference position is q:=05''35™14M255 and (5=-05°22'36".475 (J2000). 



Table 2. Results of the proper motion measurements with VERA 
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Fig. 1. — Observed spectra of the H2O maser in Orion KL. The dashed lines indicate the peak 
LSR velocity channels at the burst phase, 6.95 km s~^ and 7.58 km s~^. (a) Scalar-averaged 
cross power spectra at the burst phase on 2011 March 09 (upper) and at the quiescent phase 
on 2006 July 15 (lower). The spectra are averaged over all of the time range and all of the 
baselines, (b) Total-power spectrum (i.e. spectrum observed with each single-dish antenna 
average over the whole of time range and all of the four antennas of VERA) and scalar- 
averaged cross-power spectrum at the bursting phase on 2011 March 09. (c) Same as (b) for 
the second epoch on 2011 May 01. (d) Same as (b) for the third epoch on 2011 May 17. 
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Fig. 2. — Phase-referenced images of the bursting maser features. The delay tracking center 
position (0, 0) is a=05^35'°14M255 and (5=-05°22'36".475 (J2000). Results of the Gaussian 
fitting (major axis, minor axis, and position angle) are indicated with crosses. The syn- 
thesized beam pattern is indicated by grey ellipse at the bottom left corner in each panel, 
(a) The images at the first epoch of the bursting phase on 2011 March 09. Red and blue 
contours represent the 6.95 km s~^ and 7.58 km s~^ features, respectively. The contour levels 
are -1600, 1600, 3200, 6400, and 12800 Jy beam~^. (b) Same as (a) but for the second epoch 
on 2011 May 01. (c) Same as (a) but for the third epoch on 2011 May 17. 
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Fig. 3. — Result of position measurements with VERA. Filled circles represent the observed 
positions of the maser features. Best-fit models wit h linear proper rn otions (as indicated by 
dashed lines with arrows) and parallax of 2.39 mas (IKim et al.ll2008l ) are shown by the solid 
curves. Note that the formal errors given in Tabled] are as small as 0.03 mas, which cannot 
be shown in the figure. The reference position is the same as that in Figure [21 (a) Position 
offset in right ascension as a function of time, (b) Same as (a) for declination. Declination 
offset of feature 1 is shifted by -10 mas. (c) Movement of the bursting feature 1 on the sky. 
Crosses represent the position and the FWHM size of the maser features listed in Table [1] 
(d) Same as (c) but for the features 2 and 3. 



